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Abstract. We experimentally investigate the collapse dynamics of dipolar Bose- 
Einstein condensates of chromium atoms in different harmonic trap geometries, 
from prolate to oblate. The evolutions of the condensates in the unstable regime 
are compared to three-dimensional simulations of the Gross-Pitaevskii equation 
including three-body losses. In order to probe the phase coherence of collapsed 
condensates, we induce the collapse in several condensates simultaneously and let 
them interfere. 



PACS numbers: 03.75.Kk, 03.75.Lm 



1. Introduction 

A collapse is a fast, collective phenomenon consisting in the destruction of a multi- 
particle system happening abruptly on the time scale which governs the "usual 
dynamics". One example is the gravitational core collapse initiating a supernova. 
Happening within milliseconds, its duration is negligible compared to any of the time 
scales related to the preceding fusion stages pQ. 

In contrary to a supernova, where the experimenter is condemned to be an 
observer only, Bose-Einstein condensates (BECs) are excellent adjustable systems. 
Tailoring both the external confining potential and the interaction between the atoms 
allows to control the properties of the condensate. It is not only an ideal system to 
study questions from condensed matter physics [21 [4l [5], [6l [7l [HI [9] , but the dynamics 
of a collapse as well. 

Collapsing condensates were first observed in 7 Li [10] and 85 Rb [IT]. Both systems 
are characterized by the contact interaction, which is described by a single parameter, 
the s-wave scattering length a. From a simple model [12] minimizing the Gross- 
Pitaevskii energy functional using a gaussian ansatz for the wave function, one can 
understand the instability threshold. The energy functional consists of three terms: (i) 
the kinetic energy, which, in the sub-micro Kelvin range of almost pure condensates, 
is equal to the quantum pressure (arising from the Heisenberg uncertainty principle 
in the trap potential), (ii) the harmonic trapping potential, and (Hi) the interaction 
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energy. For repulsive contact interaction (a > 0) a global minimum exists, and the 
condensate is stable. For small enough attractive interactions (a < 0), only a local 
minimum exists at finite size; the condensate is met ast able. Finally, for sufficiently 
attractive interaction the local minimum vanishes. The potential energy of the BEC 
can be lowered without bound by contraction. Thus, the condensate is unstable. 

As the density increases, this model is too simple to describe the physics of 
the unstable cloud - atom losses due to three-body collisions have to be included. 
Being negligible at low densities, the three-body collision rate rapidly increases as the 
cloud shrinks. A three-body collision allows for the production of a dimer, where the 
third contributing atom is needed to fulfill energy and momentum conservation. The 
binding energy which is absorbed by the atom and the molecule in form of kinetic 
energy is sufficient for both to escape from the trap. Hence, instead of reaching a fully 
contracted "point-like" state, more and more atoms are lost, so that eventually the 
quantum pressure dominates over the remaining interaction energy. The dynamics 
inverts, the atoms accelerate outwards to the new equilibrium state. Because the 
three-body losses have changed the total energy, this new equilibrium state differs 
from the initial one. 

Inducing the collapse in a condensate with non-negligible dipolar interactions 
changes the dynamics completely. While the contact interaction is short-range 
and spherically symmetric, the dipolar interaction is long-range and cylindrically 
symmetric (partially attractive and partially repulsive). Therefore, dipolar 
condensates exhibit many novel phenomena even in non-collapsing systems [T3l [14j 

CEHJ CLU El CEE1 CEHJ EH EH E2JE2IESE5]. 

The anisotropy of the dipolar interaction can also be used to change the 
mechanism for the collapse - either driving or stabilizing the collapse: In a cigar- 
shaped trap with the dipoles oriented along the long axis the atoms experience 
predominantly the attractive part of the dipolar mean-field potential [26]. Hence, the 
collapse is driven by the dipolar interaction, while the contact interaction stabilizes 
the condensate. The instability occurs at positive scattering length. This situation 
is reversed if the dipoles are confined in a pancake-shaped trap. Now their dipolar 
interaction is essentially repulsive. The dipolar interaction stabilizes the condensate 
and the condensate is able to withstand negative scattering lengths for which a purely 
contact interacting BEC would have become unstable already [26J. The instability 
occurs at negative scattering length. This simple consideration strongly suggests that 
different trap geometries result in different collapse dynamics. 

In this paper we study the collapse of a chromium condensate under the influence 
of magnetic dipole interaction [27, 28, 26J. After describing our experimental procedure 
to induce the collapse, we compare the experimental data for different trapping 
potentials with 3D simulations of the Gross-Pitaevskii equation (GPE) including three- 
body losses. We show that this model provides a good description of the experimental 
data. This is not clear a priori as one might expect that the collapse dynamics 
induces many-body quantum correlations, while the GPE is a mean-field description, 
not taking correlations into account. Finally, we prove that the collapsed atomic cloud 
contains a remnant condensate by probing its phase coherence. Although simulations 
for purely contact interacting BECs have shown [29] that the observed bursts and jets 
of a collapsing condensate should be coherent, this had not been tested experimentally 
to date. 
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Figure 1. (a) Experimental setup: We realize different confining potentials by 
superimposing a one-dimensional optical lattice (green) and a crossed dipole trap 
(red). The offset coils produce the magnetic field for the Feshbach resonance, (b) 
Eddy currents retard the magnetic field B at the position of the atoms (green) 
with respect to the ramps of the offset coils (red). The time thold is defined as the 
additional holding time of the atoms in the optical trap after finishing the second 
magnetic field ramp, before the time-of-flight. (c) The corresponding scattering 
length in units of a^ g . 



2. Experimental sequence to produce collapsing dipolar condensates 

In order to create a condensate dominated by the dipolar interaction, we exploit the 
broadest of the observed Feshbach resonances in 52 Cr [30]. In the vicinity of this 
Feshbach resonance, the scattering length varies with the applied magnetic field B 
as [31] 

«{B) = a hg (l - ^^-) (1) 

where Bo w 589.1 G is the position of the Feshbach resonance, AB w (1.4 ± 0.1) G 
is its width and ab g ~ lOOaeohr is the background scattering length, with aeohr the 
Bohr radius. The magnetic field is directed along the ^-direction and determines 
the orientation of the dipoles. Using the experimental procedure described in [26], a 
condensate of approximately 50, 000 atoms is formed in a crossed far-detuned optical 
dipole trap at a ev ap ~ 85 aBohr above the Feshbach resonance. We then shape the 
external confining potentials to obtain the desired ratio of the trapping frequencies 
A := uo z /uOy by adjusting the power in the crossed dipole trap and superimposing (only 
for pancake-shaped traps A > 1) an additional one-dimensional optical lattice along the 
^-direction. The superimposed lattice is formed by two far-detuned beams (wavelength 
Aiaser = 1064 nm) crossing under a small angle ~ 8°), as shown in figure [TJ This 
results in a relatively large spacing of di at = Ai aser /[2 sin(?9/2)] = (7.4 ± 0.2) jiva 
between neighboring lattice sites. 

We first adiabatically ramp the current in the offset coils linearly in 8 ms to a 
scattering length a\ close to the point where the collapse occurs, and wait for 1 ms 
for the eddy currents (mainly due to copper gaskets in our experimental chamber) to 



special issue on cold and ultracold molecules; editors: Jun Ye and Lincoln Carr 4 



faint out. Thus, the magnetic field B(t) at the position of the atoms does not follow 
the ramp instantaneously, but must be calculated according to the equation 

B(t)=B(t)-r B B(t) (2) 

where B(t) is the magnetic field produced by the offset coils (changed linearly during 
ramps) and t# w (0.57 ± 0.05) ms is the 1/e lifetime of the eddy currents, measured 
by Zeeman spectroscopy [32]. After 1 ms of waiting time at ai, we start a second ramp 
from cl[ to <2f < a cr it(A), where a cr it(A) is the critical scattering length for the given 
trapping potential [26], so that the collapse occurs. We hold the atoms in the trap for 
an additional time thoid a t <2f before releasing them and taking a time-of-flight image. 
In order to get the maximal absorption cross section, we split the time-of-flight into 
two parts: a first one, lasting 4 ms, at the magnetic field corresponding to <2f (in order 
not to disturb the dynamics) and a second part, lasting again 4 ms, where the large 
magnetic field along z is replaced by a field of 11 G along the x-direction. We checked 
that this procedure does not disturb the image. 

The observed integrated density distribution is bimodal and consists of a broad 
isotropic thermal cloud and an anisotropic remnant BEC (see e.g. Fig 1(b) of [32J). 
Because the size of the thermal cloud as well as its atom number does not depend on 
^hoid, it is unlikely to contribute to the collapse dynamics. In the following, we have 
subtracted it from the images to increase the contrast. 

3. Collapse dynamics for different trap geometries 

The anisotropic character of the dipolar interaction (dipoles side-by-side repel each 
other, while dipoles in a head-to-tail configuration attract each other) has a strong 
effect on the stability of a dipolar condensate: varying the geometry of the confining 
harmonic trap from prolate to oblate (the symmetry axis being the one along which 
the dipoles are aligned) stabilizes the condensate, as was demonstrated experimentally 
in [26]. Our previous experimental study of the collapse dynamics of a dipolar 
condensate [32] was restricted to an almost spherical trap; in the following we study 
the influence of the trap geometry on the collapse dynamics. 

3.1. Prolate traps 

For purely contact interacting condensates the time scale which governs the "usual 
dynamics" is set by the largest trap frequency. In contrast, for dipolar condensates 
this time scale is given by the largest radial trap frequency r := 1 //radial, because the 
collapse is induced in this direction [32] . 

Figure [2] presents the collapse of a dipolar condensate in a very elongated prolate 
trap (A ~ 0.12). Here we ramp from an initial scattering length a\ — (35 ± 2) aBohr to 
af = (8 ± 3)aBohr, which lies below the critical scattering length a crrt ~ 12aBohr- On 
the time scale of the radial trap period r « 0.76 ms the condensate only starts to split. 
This is explicitly shown in figure [2] (d): The atom number does not drop "abruptly" 
to its final value, but instead it changes linearly on the time scale r. Therefore, in the 
case of figure [2] we observe only a "moderate" collapse. 

Each panel in figure [2] (a) and (b) presents the integrated column density 
J dx\i(j(r,t)\ 2 of an absorption image after 8 ms of free expansion. While the upper row 
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Figure 2. Comparison of the experimental absorption images (a) with the 
simulations (b) for a cigar-shaped trap with trapping frequencies (fx, fy, fz) = 
(1312, 1311, 161) Hz and trap ratio A « 0.12. Each image shows the averaged 
column density of five pictures after 8 ms of time-of- flight. The field of view is 
250 X 250 jum 2 . The final value of <2f was (8 ± 3) aBohr (also in the simulation). 

(c) Simulated in-trap absorption images. The field of view is (y, z) = (4.8, 25) jum. 

(d) The remnant condensate atom number for different holding times in units of 
the radial trap period r = l// ra diai — 0-^6 ms. The solid line is the result of the 
numerical simulation (see text), without any adjustable parameter. 



shows the experimental data, the lower rows is obtained from a numerical simulation 
of the three dimensional Gross-Pitaevskii equation 

ih^(r,t) = J - ^-V 2 + V tTap + f U(f- f ,t)|V(r',t)| 2 dV - ^\mt)A W,t) 
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Figure 3. (a, b) Collapse dynamics for different holding times thold in a weakly 
cigar-shaped trap with trap frequencies (fx, fy, fz) ~ (650, 520, 400) Hz, trap ratio 
A « 0.7 and radial trap period r ~ 1.7 ms. Each picture is a single absorption 
image after 8 ms of time-of-flight and has its own optical density scale. The field 
of view is 250 x 250 /im 2 . (c) Simulated in-trap absorption images. The field of 
view is 6.9 x 6.9/xm 2 . (d) N-qbc versus thold- The solid line is the simulation 
result. 



with the contact and dipolar interactions 



U(f, t) = 4nn l a{t) S(r) + ^! 1 - 3 7 2 



47T |f] 2 

where m is the atomic mass, /iq the magnetic permeability of vacuum, \i = 6/iBohr the 
magnetic moment of chromium, and is the angle between r and the magnetic field 
B || e z . The imaginary term describes the atom losses. The simulations contain no 
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free parameter, as the three-body coefficient L3 = 2 x 10 -40 m 6 /s was estimated from 
measurements. Details about the simulations are described in [32] . 

Figure [3] shows the collapse in a cigar-shaped trap with trapping frequencies 
(fx, fy, fz) ~ (650, 520, 400) Hz, corresponding to A ^ 0.7, and radial trapping period 
r ~ 1.7 ms. We start with TVbec = 13,500 ± 1,500 atoms before ramping the 
scattering length from a\ = (35 ± 2) aBohr to af = (8±3)aBohr, which lies 4aBohr 
below the critical scattering length. The absorption images indicate three different 
stages: First, for thoid = ms, the condensate is strongly elongated along the magnetic 
field direction, demonstrating strong dipolar interactions [25] . Second, we observe a 
change of ellipticity after 0.4 ms ~ 0.3 r. This is a consequence of the radial implosion 
and subsequent explosion (a stable, cigar-shaped dipolar condensate does not invert 
its ellipticity during the free expansion [28]). Third, we observe a splitting of the 
condensate in axial direction as in figure [21 but now after 0.8 ms « 0.5 r. For longer 
holding times the splitting becomes more prominent, but the dynamics is already 
completed after half a trapping period. This is shown in figure [3] (d). The remnant 
atom number A^bec(^) drops from 12,000 to 4,000 within this time scale. Note that, 
for this figure (and only this one), the simulation results agree better with the data 
when performed with af = 2 aeohr (the results shown in figure [3] are for this value of 
af). This slight discrepancy is most probably due to a slow drift of the magnetic field 
that occurred between the data acquisition and the calibration of the scattering length 
(the calibration procedure, described in details in [25] . requires the accumulation of a 
large quantity of data). 

3.2. Oblate trap 

As discussed, different trap geometries are expected to result in different collapse 
dynamics. So figures [2] and [3] have to be compared to the time evolution of an almost 
spherical and a pancake-shaped trap. While the collapse in an almost spherical trap 
is published in [32] , we present the pancake-shaped trap in figure [U 

The pancake-shaped trap is formed by superimposing two additional lattice beams 
onto a condensate, which was produced in the crossed dipole trap, see figure [TJ 
Depending on the non-stabilized relative phase of the two lattice beams, it is in 
principle possible that the condensate splits into two. Practically, we have never 
observed interference fringes, even at expansion times long enough, so that the fringe 
spacing ht to f / (md\ a t) was larger than our 6 fim resolution limit. Switching off the 
trapping potential immediately after finishing the magnetic field ramp results in a 
very elongated condensate, in agreement with the simulations (figured]). 

The cylindrical symmetry of the pancake-shaped trap allows us to recover the 
three-dimensional density distribution n(p, z) from the two-dimensional absorption 
image n a b s (2/, z). Figured] (d) and (e) show isodensity surfaces for two different 
densities obtained from the inverse Abel transformation [33] 



for thoid = 0.4 ms, where Jo is the Bessel function of the first kind. As in the 
almost spherical trap [32J, the collapsed cloud exhibits a d-wave symmetry. While 
the isodensity surface for "high" densities contains a central disc and two separated 
"blobs" along the symmetry axis, the three parts merge for "low" densities. Because 
the Abel transformation is very sensitive to noise, we cannot extract in a reliable way 
the kinetic energy released during the collapse from our images. 
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Figure 4. (a, b) The collapse dynamics for different holding times thold in 
a pancake-shaped trap, which corresponds to trap- frequencies (fx, fy, fz) ~ 
(400,400,3400) Hz, trap ratio A w 8.5 and radial period r ~ 2.5 ms. The upper 
row present the average of five absorption images after 8 ms of free expansion. 
The final ramp starts at a\ = (30 zb 2) a-Bohr an d stops at a,f = ( — 13 ± 2) a-Bohr- 
The simulations give a cr it = ( — 1.5 zb 0.5) a-Bohr- ( c ) Simulated in-trap absorption 
images. The field of view is (y, z) = (8.9,4.5) fim. Isodensity surfaces for "high" 
(d) and "low" (e) densities of the thold — 0-4 ms image obtained from the Abel 
transformation. 



4. Testing the coherence of the collapsed cloud 

By analyzing the collapse of condensates confined in different trap geometries, we 
showed that the survival density pattern exhibits two parts: one part, which is well 
described by a thermal cloud and a second part, which, due to its high optical density, 
we interpreted as a remnant condensate [32]. In order to confirm this interpretation, 
we checked the coherence of interfering condensates in the case of pancake-shaped 
traps. 

For that purpose, we produce several condensates by superimposing the two 
beams for the optical lattice before finishing the evaporation to reach quantum 
degeneracy. Here the extension of the cloud is still large enough, so that the atoms 
occupy three to five adjacent lattice sites. The exact number of occupied sites changes 
from shot to shot, depending on the relative phase of the two lattice beams, which 
are not actively stabilized. After condensing, the BECs on different lattice sites have 
random phases with respect to each another. The single particle tunneling rate is 
vanishingly small (< 10 _30 /s), so that no phase coherence between adjacent sites is 
built up on the time scale of the experiment. 

It is known [34] that the interference fringes are not washed out if few independent 
condensates interfere instead of only two. Taking the same trap geometry as in figure[H 
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Figure 5. Interference pattern of three to five independent condensates 
for different holding times thold m the pancake-shaped trap, (fx,fy,fz) ~ 
(400, 400, 3400) Hz. The lower graphs show the y- integrated column densities 
Y^ y J dx\i\)\ 2 (time-of-flight 18 ms, field of view 110 x 690 /im 2 ). 

but extending the time-of-flight to 18 ms, we obtain the interference patterns shown 
in figure [U As expected, the absolute position of the interference fringes changes from 
shot to shot, but they can be clearly seen on each image. 

While for holding times shorter than 0.2 ms or longer than 0.5 ms interference 
fringes with a high contrast are visible, we observe no interference fringes at the center 
of the two clouds for 0.2 ms < thold < 0.5 ms. A possible interpretation is as follows: 
For thold = ms no collapse occurs and the observed fringes are similar to those 
of two point sources. For thold = 0.2 ms and 0.4 ms the condensates do collapse, 
but this happens during the time-of-flight and after the clouds overlapped: e.g. for 
^hoid = 0.4 ms the condensates start to overlap for t to f = 0.4 ms, but the collapse 
happens at t to f = 0.8 ms. Probably this induces a complicated phase distribution in 
each of the condensate and integration over the line of sight washes out the interference 
fringes. On the other hand, if the collapse happens in-trap (e.g. for thold = 0.8 ms ) 
the fringes are formed by the remnant condensates. Again we recover the fringes as if 
the atoms would belong to two coherent point sources. 

5. Conclusions 

We have investigated theoretically and experimentally the collapse dynamics of dipolar 
condensates in prolate and oblate harmonic trapping potentials. As expected, the 
collapse dynamics depends on the trap geometry, although the qualitative behavior 
is similar for different traps. The simulations containing no adjustable parameter 
reproduce the experimental results well. By simultaneously inducing a collapse in 
several condensates and let them interfere, we showed that the collapsed cloud contains 
a coherent remnant condensate. 

A clear direction of further studies is to use the interferometric technique of 
section |4] to obtain experimental evidence for the vortex rings predicted in [32J . The 
contrast of the interference fringes is not high enough to support their evidence yet. 
An other interesting extension is the study of two-dimensional solitons [35, 36J, which 
are expected to appear just above the instability threshold. 
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